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SECTION  I 


INTRODUCTION 

Skin  friction  can  be  significantly  larger  for  turbulent  flow 
over  a  rough  surface  as  compared  with  an  equivalent  flow  over  a 
smooth  surface.  Many  systems  of  engineering  interest,  such  as 
re-entry  vehicles,  heat  exchangers,  aircraft,  turoines  and  piping 
networks,  have  surfaces  which  are  often  rough  in  the  aerodynamic  or 
hydraulic  sense.  Therefore,  there  is  significant  interest  in 
accurate  predictive  models  for  turbulent  flows  over  rough  surfaces. 

In  turbulent  flow  analysis,  use  of  time-averaged  equations  leads 
to  the  necessity  of  formulating  a  turbulence  model  with  empirical 
input  to  achieve  closure.  A  similar  situation  exist  in  the  analysis 
of  flow  over  rough  surfaces.  Unless  the  equations  can  be  solved  in  a 
grid  that  is  fine  enough  to  resolve  the  surface  roughness  geometry,  a 
roughness  model  with  empirical  input  is  necessary. 

When  considering  the  development  of  predictive  models  for 
turbulent  flow  over  rough  surfaces,  it  is  necessary  to  divide  the 
surfaces  into  two  broad  classes  based  on  the  geometry  of  the 
roughness  elements.  These  two  classes  are  three-dimensional 
distributed  roughness  and  two-dimensional  rib-type  roughness.  The 
distributed  roughness  is,  as  its  name  implies,  composed  of  a 
collection  of  small,  mostly  three-dimensional  protrusions  such  as  an 
array  of  hemispheres  or  the  texture  of  a  casting.  Likewise,  rib-type 
roughness  is  composed  of  two-dimensional  strips,  such  as  in  a  fiber- 
reinforced  composit.e.  While  related,  the  interactions  of  these  two 
classes  of  roughness  with  the  flow  field  involve  different  physical 
behaviors  and  require  separate  treatments. 

This  report  presents  a  discrete  element  model  for  turbulent  flow 
over  rib-type  roughness.  A  brief  discussion  of  the  background  of 
roughness  modeling  is  given  below  'nd  is  followed  by  a  description  of 
the  current  model  and  comparisons  of  this  model  with  experimental 
data  sets.  A  more-  detailed  discussion  of  roughness  modeling  can  be 


SECTION  II 


BACKGROUND  OF  ROUGHNESS  MODELING 


The  first  in-depth  treatments  of  the  effects  of  surface 
roughness  on  fluid  flow  were  those  of  Nikuradse  [2]  and  Schlichting 
[3].  Schlichting  proposed  that  rough  surfaces  be  classified  based  on 
|  equivalent  sand-grain  roughness  (k3).  He  defined  ks  as  the  size  of 

sand  grain  in  Nikuradse's  pipe  flow  experiment  which  would  give  the 
same  skin  friction  as  that  observed  on  a  particular  rough  surface, 
and  he  experimentally  determined  ks  for  a  variety  of  surfaces. 
Schlichting's  work  was  for  many  years  the  accepted  standard,  and  much 
of  the  subsequent  theoretical  and  experimental  research  incorporated 
his  results,  explicitly  or  implicitly,  into  models  and  phenomeno¬ 
logical  explanations.  Because  Schlichting's  data  reduction  method 
|  was  flawed  (mainly  by  neglecting  the  effects  of  the  side  walls  on  his 

fi:  1  aspect  ratio  channel),  much  of  the  subsequent  work  which 
referenced  his  results  and  attempted  to  cast  various  roughnesses  in 
terms  of  sand-grain  equivalent  are  also  in  error.  Coleman  et  al.  [Ji] 
|  reported  the  flaw  in  Schl ichting's  data  reduction  procedure  and 

presented  corrected  values  for  his  data. 

However,  even  with  a  consistent  data  base,  problems  remain  when 
attempts  are  made  to  extend  the  sand-grain  roughness  concept  to 
specific  roughness  geometries.  Implicit  in  the  proper  use  of 
equivalent,  sand-grain  roughness  are  the  dual  assumptions:  (1) 
roughness  effects  can  be  adequately  described  by  a  single  parameter 
and  (2)  over  a  wide  range  of  conditions  the  roughness  behaves  as  if 
it  were  sand  grain.  As  Reynolds  [b]  points  out  these  are  not  viable 
assumptions  for  many  roughness  configurations.  Another  problem  with 
this  approach  is  how  to  determine  k3  when  no  skin  friction  data  exist 
for  a  particular  surface. 

The  typical  tactic  taken  is  to  try  to  relate  ks  to  geometrical 
attributes  of  a  specific  roughness  through  a  correlation.  Dvorak  [6], 
Simpson  [7]  and  Dirling  [8]  have  presented  ks/ka  as  a  function  of 
various  geometrical  parameters.  In  these  correlations,  ks  is  the 
equivalent  sand-grain  roughness ,  and  k,  is  the  nominal  roughness 
height .  The., 1  correlations  do  not  correlate  the  data  well,  and  they 


rely  primarily  on  Schlichting's  results  which  the  authors  have  shown 
to  bo  in  error  ['0.  Figure  1  presents  the  correlation  of  Dirling 
along  with  th^  original  data  ho  used.  In  this  figure,  A  is  the 
roughness  density  parameter,  As  and  Ap  are  the  windward  roughness 
element  surface  area  and  the  projected  windward  surface  area 
respectively,  and  lr  is  the  average  center-to-center  spacing.  While 
trie  Pirling  correlation  tends  to  yield  the  correct  trends,  the 
scatter  is  quite  severe  as  ks/ka  varies  by  a  factor  of  three-to-four 
for  most  A  values.  Berg  [9],  foe  example,  found  that  the  Dirling 
correlation  gave  a  ks/ka  value  three  times  larger  than  his  experimen¬ 
tally  determined  values. 

The  correlation  also  treats  two-dimensional  and  three-dimen¬ 
sional  roughness  in  the  same  manner,  even  though  the  nature  of  the 
flow  around  the  two  types  is  very  different.  This  is  evident  also 
from  Figure  1 ,  where,  except  for  the  most  densely  packed  spherical 
roughness  elements  of  Schl ioht.ing,  the  two-  and  three-dimensional 
roughness  correlations  are  distinct.  The  flow  over  two-dimensional 
elements  is  primarily  dependent  on  the  spacing/height  ratio  whereas 
spacing,  shapes  and  height  are  all  important  for  flow  over  and 
around  three-dimensional  elements.  Perry  et  al .  [10]  found  that  for 
transverse  bar  elements  with  certain  height-to-spacing  ratios  the 
skin  friction  does  not  scale  with  the  roughness  height.  The  ks 
concept  is  particularly  poor  for  rib-type  roughness.  Morris  [11] 
presented  a  classification  of  rib-type  roughness  which  was  based 
primarily  on  t  lie  rib  spacing  to  height  ratio  (L/k) .  Reynolds  [5] 
identifier.  t  >ur  distinct  relationships  of  skin  friction  with  Reynolds 
number ,  ..nly  tir  of  which  is  "sand  grain"  in  behavior.  Several 
rre’a'i-  i.::  !  r  fi  i<  M  .n  coefficient  have  been  presented  based  on 
tic.  •  it-':  ,n:.t  •  i  •  r  ••  pjivalent  sand-grain  roughness  (Han  et  al. 

.•••  ri olive  nature  of  the  implications  of  basing 
a  r  u.n. i,-.  ‘  r.  :■  .  •  i  1  .r--ly  on  sand  grain  or  equivalent  sand-grain 

i  -  igtiii-' . :  *  r.<-  mom  funriunent.a  1  ly  based  discrete  element 

rougt.n-  :c.  Ci  )•. .  w  c>  ii  car..  This  .ipproaoh  models  the  interaction  of 
t  r,<  roughness  n.d  tn*  1  l.>w  by  account  i  r.g  for  the  blockage  effect  of 
tti<  --oughness  -Hi'!  *  he  drag  force  which  the  roughness  elements  exert 


Equivalent 


on  the  fluid.  Roberson  and  his  co-workers  (Reference  13.  for 
example)  used  some  aspects  of  this  model  in  their  insightful  work. 
Finson  [  1  -I ] ,  L.in  and  Bywater  [15]  and  Christoph  and  Pletcher  [16] 
have  reported  work  whi'-h  uses  these  concepts.  The  most  recent 
results  have  been  presented  by  the  authors  [1,  17].  The  discrete 
element  roughness  model  for  rib-type  roughness  is  examined  in  the 
next  section. 
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Figure  1.  Effective  Sand-Grain  Roughness 
Dir  ling  [8] 
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SECTION  III 


DISCRETE  ELEMENT  MODEL  FOR  RIB-TYPE  ROUGHNESS 

It  has  long  been  recognized  that  the  governing  phenomenon  of  the 
interaction  between  rib  roughness  and  the  flow  is  the  momentum  sink, 
which  results  in  the  form  drag  of  the  roughness.  Morris  [11] 
suggested  that  form  drag  concepts  could  be  used  to  predict  the 
friction  factor  for  flow  over  rib-type  roughness.  Lewis  [18]  used 
these  ideas  as  a  tool  in  the  interpretation  of  his  experimental  data. 
However,  a  complete  discrete  element  formulation  has  not  been  used  to 
predict  the  turbulent  flow  field  over  surfaces  with  rib-type 
roughness . 

The  discrete  element  approach  considers  the  mass  and  momentum 
transport  processes  on  the  collection  of  individual  roughness 
elements  and  the  smooth  surface  between  them.  The  basic  idea  of  the 
discrete  element  approach  is  to  formulate  a  system  of  partial 
differential  equations  which  describe  the  mass  and  momentum  transport 
over  and  between  the  elements.  In  this  method  the  roughness 
effects  (form  drag  and  blockage)  are  taken  as  an  integral  part  of  the 
problem. 

The  development  and  validation  of  the  discrete  element  roughness 
model  of  Coleman  et  al.  is  discussed  in  detail  in  References  1  and 
17.  An  abbreviated  discussion  is  presented  herein.  In  the  discrete 
element  approach  the  rough  surface  is  assumed  to  be  composed  of 
distinct  roughness  elements.  This  model  is  based  on  consideration  of 
the  physics  of  the  interactions  between  the  discrete  roughness 
elements  and  the  flow  and  does  not  depend  on  the  assumption  of  any 
sand-grain  roughness  equivalent.  An  array  of  two-dimensional 
roughness  elements  perpend icular  to  the  streamwise  direction,  x,  is 
used  in  the  following  discussion. 

The  differential  equations  including  roughness  effects  are 
derived  by  applying  the  basic  conservation  statements  for  mass  and 
momentum  to  a  control  volume  (CV)  such  as  that  shown  in  Figure  2. 

This  CV  is  shown  with  an  exaggerated  length,  6x,  in  the  primary  flow 


Figure  2.  Control  Volume  for  Flow  Over  Two-Dimensional  Roughness 

penetrated  by  roughness  elements.  The  physical  effects  of  the 
roughness  elements  on  the  fluid  in  the  CV  are  modeled  by  considering 
the  flow  blockage  and  by  postulating  that  the  total  force  of  the 
elements  on  the  flow  can  be  incorporated  as  a  drag  force. 

Basic  to  this  approach  is  the  idea  that  the  two-dimensional, 
time-averaged  turbulent  boundary  layer  equations  can  be  applied  in 
the  flow  region  below  the  crests  of  the  roughness  elements. 
Therefore,  the  flow  variables  must  be  viewed  as  having  been  averaged 
over  the  transverse  (z)  direction  and  averaged  in  the  longitudinal 
(x)  direction  over  an  appropriate  x  distance. 

Mass  and  momentum  balances  were  made  for  this  control  volume  in 
order  to  derive  the  partial  differential  equations  that  describe  the 
flow  over  a  rough  surface.  The  resulting  partial  differential 
equations  were  then  time  averaged  (Reynolds  averaged)  and  the  bound- 
ary-layer  assumptions  invoked.  The  resulting  equations  are 
Cont  i  nu  i  t.  y 


x-Momentum 


where  Bx  and  By  are  the  blockage  factors,  L  is  the  distance  between 
ribs  and  Cp  is  the  local  element  drag  coefficient. 

A  comparison  of  the  above  equations  with  the  standard 
time-averaged  turbulent  boundary-layer  equations  for  smooth  walls 
reveals  the  results  of  including  the  discrete  roughness  elements.  A 
term-by-term  examination  of  the  differences  is  appropriate.  The 
parameter  B  represents  the  fraction  of  the  control  surface  (CS)  open 
to  the  flow.  Bx  is  the  fraction  of  the  CS  perpendicular  to  the 
x-coordinate  that  is  open  to  the  flow;  By  is  the  fraction  of  the  CS 
perpendicular  to  the  y-coordinate  that  is  open  to  the  flow.  Bx  and 
By,  the  blockage  factors,  also  appear  in  the  x-momentum  equation. 

The  additional  term  in  the  x-momentum  equation  represents  the  form 
drag  of  the  roughness  elements  and  is  represented  in  terms  of  a  drag 
coefficient.  This  term  acts  as  a  distributed  momentum  sink. 

The  skin  friction  coefficient  for  the  discrete  roughness  element 
model  is  as  follows: 


Cf  = 


(6x)w  TW  +  -  —  „  pCDU2(dy) 


(5) 


The  subscript  w  indicates  values  evaluated  at  the  wall,  and  the  e 
indicates  boundary-layer  edge  conditions.  The  smooth  wall  shear 
stress  tw  is  given  by  u(3U/9y)w.  The  skin  friction  coefficient  is  a 
non-dimensional ized  form  of  the  apparent  wall  shear  stress.  The 


apparent  wall  shear  stress  is  composed  of  the  viscous  shear  term 
acting  over  the  portion  of  the  wall  not  occupied  by  roughness 
elements  (the  smooth  portion)  plus  the  resultant  form  drag  on  the  rib 
roughness  elements.  The  proper  formulations  for  the  blockage 
factors,  8X  and  By,  and  the  drag  coefficient,  Cp,  along  with 
turbulent  closure,  complete  the  model. 


ROUGHNESS  MODEL  FOR  RIB-TYPE  ROUGHNESS 


For  three-dimensional  roughness,  the  blockage  factors  required 
for  the  discrete  element  model  of  Coleman  et  al.  [1,  17]  were 
functions  of  only  the  roughness  element  geometry.  However,  the 
blockage  factors  for  the  discrete  element  formulation  for 
two-dimensional,  rib-type  roughness  require  some  fluid  mechanical 
considerations. 

Figure  3,  taken  from  Lewis  [18],  provides  a  phenomenological 
schematic  of  the  salient  features  of  the  interaction  of  rib-type 
roughness  and  the  turbulent  boundary  layer  in  which  it  is  immersed. 


FULLY  TURBULENT 

FLOW  CORE 

- » 


HP  =  High  Pressure,  LP  =  Low  Prossuro 


Figure  3.  Schematic  of  Salient  Features  of  Flow  Over  Two-Dimensional 
Roughness  (From  Lewis  [18]) 

The  dominant  features  of  the  interaction  are  flow  separation  and 
reattachment.  These  regions  of  separated  flow  give  rise  to  a 
momentum  sink  which  results  in  the  form  drag  of  the  element.  Three 
separated  regions  are  present  on  a  given  rib:  (1)  a  large  (in 
comparison  to  the  rib  size)  region  downstream  of  the  element,  (2)  a 
smaller  region  over  the  upstream  face  of  the  element,  and  (3)  a  very 


small  (in  comparison  to  the  rib  size)  region  on  the  upstream  side  of 
the  top  of  the  element.  The  windward  (upstream)  face  sees  a  much 
higher  pressure  level  than  the  leeward  (downstream)  face. 

If  the  rib  spacing  is  large  enough,  the  boundary  layer 
reattaches  to  the  wall  at  some  point  between  the  ribs.  The 
quantitative  aspects  of  these  separations  and  reattachments  are 
highly  dependent  upon  the  rib  geometry  and  attributes  of  the 
turbulent  boundary  layer.  In  situations  where  the  streamwise  flow 
reattaches  to  the  wall,  the  streamlines  drop  below  the  crests  of  the 
two-dimensional  ribs  and  then  move  above  the  crests  in  the  neighbor¬ 
hood  of  the  following  elements. 

Quantitative  details  in  the  form  of  static  pressure  isobars  for 
the  flow  field  over  a  square-rib  configuration  are  presented  in 
Figure  4.  This  figure  was  taken  from  Lawn  [19]  and  is  illustrative 
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Figure  h .  Isobars  of  Typical  Flow  Over  Two-Dimensional  Rib  Roughness 
(from  Lawn  [19]) 

of  the  comp'  xity  of  the  flow  structure  encountered  in  turbulent 
boundary  layer  flow  over  two-dimensional  roughness  configurations. 

The  pressure  variation  over  the  upstream  face  of  the  element  shows  a 
decrease  from  the  base  of  the  element  to  the  top  of  the  element; 
whereas  the  pressure  on  the  rear  face  (downstream  side)  of  the 
element  is  at  a  much  lower  level  and  is  essentially  constant. 


The  information  presented  in  Figures  3  and  4  is  indicative  of 
the  complexity  of  the  viscous  flow  over  two-dimensional  roughness. 

The  presence  of  the  boundary  layer  is  responsible  for  much  of  the 
complexity.  Good  and  Joubert  [20]  point  out  that  while  bluff  bodies 
immersed  in  an  external  flow  have  essentially  a  constant  form-drag 
coefficient  (for  all  Reynolds  numbers  above  a  few  thousand),  bluff 
bodies  immersed  in  a  turbulent  boundary  layer  have  form-drag 
coefficients  that  are  dependent  upon  the  Reynolds  number  because  of 
the  distribution  of  the  incident  momentum  flux  and  the  behavior  of 
the  separated  regions. 

If  the  flow  about  two-dimensional  rib  roughness  is  to  be 
calculated  in  detail,  then  the  Navier-Stokes  equations  with  an 
appropriate  grid  and  turbulence  model  must  be  used.  For  roughnesses 
with  small  spacings,  a  grid  system  which  resolves  the  surface 
geometry  presents  a  difficult  problem.  Indeed,  the  complexity  of  the 
physics  of  the  rib  roughness/turbulent  boundary  layer  interaction  is 
so  great  that  only  recently  have  turbulence  models  and  numerical 
techniques  been  developed  that  can  adequately  predict  important 
features  of  the  flow.  The  recent  work  of  Benodekar  et  al.  [21]  using 
a  modified  two-equation  (<,  e)  turbulence  model  is  typical  of  recent 
(1985)  predictive  efforts.  While  they  predicted  important  nuances  of 
some  characteristics  of  the  flow  field  over  a  single  rib  element, 
some  anomalies  were  present.  For  example,  the  predictions  of  the 
upstream  pressure  rise  were  in  excellent  agreement  with  the  data,  but 
the  prediction  of  the  downstream  pressure  variation  showed  some 
discrepancies.  The  predicted  velocity  profiles  were  not  in  agreement 
with  measured  profiles  downstream  of  the  roughness  element. 

The  preceding  background  information  is  sufficient  to  develop 
the  roughness  model.  Equations  (1)  and  (2)  were  obtained  by 
averaging  the  effects  of  roughness  in  the  longitudinal  (x)  and 
transverse  (z)  directions.  Thus,  the  blockage  terms  and  the  form 
drag  (momentum  sink)  expression,  as  they  appear  in  equations  (1)  and 
(2),  are  in  effect  expressed  on  a  per  unit  length  basis  in  the 
streamwise  direction.  Hence  the  discrete  element  roughness  model  for 
two-dimensional,  rib-type  roughness  makes  no  attempt  to  model  details 
of  the  flow  field  at  discrete  roughness  locations,  but  instead  models 
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the  effects  of  blockage  and  form  drag  in  such  a  manner  that  the 
average  wall  resistance  (viscous  shear  plus  form  drag)  is  accurately 
predicted.  The  regions  of  separation,  which  would  require 
Navier-Stokes  computations,  are  replaced  by  the  roughness  model  so 
that  parabolic  (marching)  computational  techniques  with  the  turbulent 
boundary  layer  equations  can  be  utilized. 

In  a  manner  similar  to  that  of  Morris  [11],  the  flow  over 
rib-type  roughness  is  divided  into  three  regimes: 

1 )  Reattached  flow  -  the  separated  region  behind  the  rib 
reattaches  to  the  smooth  surface  between  the  ribs;  viscous 
skin  friction  is  a  significant  factor. 

2)  Unreattached  flow  -  the  separated  regions  behind  and  in 

front  of  the  ribs  merge;  viscous  skin  friction  is  relatively 

unimportant . 

3;  Skimming  flow  -  the  region  between  the  ribs  is  completely 
filled  with  a  recirculating  flow,  the  surface  appears 
semi-smooth. 

The  formulation  presented  herein  applies  only  to  regimes  1  and  2. 
Figure  5a  and  5b  schematically  illustrate  the  models  for  these 
regimes.  The  separated  region  upstream  of  the  roughness  is  taken  to 
be  the  3ame  length  as  the  element  height.  This  agrees  qualitatively 
with  the  interaction  given  by  Lawn  [19]  and  presented  in  Figure  4 , 

and  with  the  suggestion  of  Lewis  [18].  For  the  purposes  of  blockage 

factor  calculations,  the  separated  region  i3  taken  to  be  that  area 
enclosed  by  a  line  extending  from  the  separation  location  to  the  tip 
of  the  element. 

The  downstream  separation  length  is  taken  as  Yk  and  will  be 
discussed  in  more  detail  later.  Similar  to  the  upstream  region,  the 
downstream  separated  region  is  taken  as  extending  from  the  tip  of  the 
roughness  element  to  the  reattachment  point.  The  factor  6  is  defined 
as  that  fraction  of  a  control  surface  that  is  open  to  flow  and  the 
action  of  viscous  shear  due  to  attached  streamwise  flow.  The  large 
separated  regions  before  and  after  the  rib  act  as  blockage  to  the 
main  flow  field.  Referring  to  Figure  5a  the  blockage  factor  for  an 


x-z  plane  is  seen  to  be 


6y  =  1  -  d/L  (6) 

where  d  =  w  +  (Y  +  1)(k  -  y) 

The  blockage  factor  for  the  y-z  plane,  Bx>  is  more  difficult  to 
define.  For  some  y-z  planes  the  flow  is  completely  blocked  below  the 
crest  of  the  ribs.  For  others  there  is  no  blockage.  If  an  average 
is  taken  over  a  length  containing  several  ribs,  an  average  blockage 
factor  in  y-z  planes  is  found  which  is  approximately  equal  to  the  x-z 
blockage  factor  given  in  equation  (6).  Above  the  ribs  the  blockage 
factor  is  1  for  both  8X  and  By 

For  flow  where  the  separated  regions  merge,  regime  2,  the  same 
equations  apply.  However,  since  the  wall  is  completely  covered  by 
the  recirculating  flow,  a  new  origin  for  y  must  be  taken.  Referring 
to  Figure  5b  and  following  Lewis  [18],  the  new  origin  is  located  a 
distance  A  above  the  base  wall  where 

A  =  k  -  (L  -  w)/(Y  +  1)  (7) 

The  effect  of  this  construction  is  to  treat  the  roughness  as  shorter 
ribs  in  regime  2  flows. 

As  shown  by  the  data  of  Lawn  [19]  and  presented  in  Figure  i»,  the 
fluid  flowing  over  the  roughness  ribs  results  in  nearly  constant 
pressure  on  the  leeward  side  of  a  rib.  The  pressure  on  the  windward 
side  decreases  from  the  base  to  the  top  of  the  rib.  It  is  this 
windward-to-leeward  side  pressure  difference  which  results  in  the 
local  form  drag.  As  shown  in  equation  (2)  this  local  drag  force  is 
ca3t  in  terms  of  a  local  drag  coefficient 

CD  -  2AP/ pU^  (8) 

."ince  the  pressure  difference  decreases  with  distance  above  the  base 
wall  while  the  velocity  increases  with  distance  above  the  base  wall, 
Cp  must  be  a  decreasing  function  of  the  distance  above  the  base 
surface.  Based  on  Lawn's  [19]  detailed  pressure  and  velocity 
distribution  measurements  in  the  neighborhood  of  roughness  ribs,  the 
following  formulation  was  used 


The  constant  Cqq  was  found  by  calibration  using  comparisons  of 
numerical  experiments  with  base  data  sets. 

For  turbulent  flow  some  closure  model  must  be  adopted  to  relate 
the  shear  stress,  t,  to  the  mean  flow  field.  Here  the  Boussinesq 
approximation  is  invoked: 

dU 

t  =  (p  +  pt)  —  (10) 

3y 

The  eddy  viscosity,  p-p,  is  modeled  using  the  usual  smooth  wall 
Prandtl  mixing-length  model  with  van  Driest  damping.  Since  the 
roughness  effects  were  incorporated  entirely  into  the  roughness 
model,  the  mixing-length  model  used  was  the  conventional  one  used  for 
smooth  surfaces.  That  model  is  (Schetz  [22]) 

im  =  0.H0y[l  -  exp(  -y  +  /A+ )  ]  ;  8,m  <  0.09  6 

Hm  =  0.096  otherwise.  (11) 

The  appropriateness  of  using  a  smooth-wall,  mixing-length  model  was 
validated  by  Taylor  et  al.  [1]  for  three-dimensional  roughness. 

The  time-averaged  turbulent  boundary-layer  equations  as 
formulated  for  discrete  element  roughness  and  with  the  aforementioned 
turbulence  and  roughness  models  were  solved  numerically  using  an 
implicit  finite  difference  technique  similar  to  that  of  Adams  and 
Martindale  [23].  Comparisons  of  the  discrete  element  roughness  model 
as  developed  herein  with  several  experimental  data  sets  are  presented 
in  the  next  sections. 
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SECTION  V 


COMPARISON  OF  PREDICTIONS  WITH  EXPERIMENTAL  DATA 

I 

Comparisons  have  been  made  between  predictions  using  the  model 
described  herein  and  several  experimental  data  sets  [12,  24,  25,  26]. 
The  range  of  Reynolds  numbers  based  on  diameter  or  hydraulic  diameter 
|  for  these  comparisons  was  approximately  6,000  to  200,000.  The  data 

show  [11,  18]  that  the  length  of  the  leeward  separated  region  is 
between  3k  and  8k.  For'  all  the  predictions  shown  in  this  paper  Y  =  5 
was  used.  This  value  gave  the  best  results  for  the  more  closely 
|  spaced  rib  elements,  L/k  <  10.  For-  the  more  widely  spaced  rib 

elements,  the  model  was  relatively  insensitive  to  Y  within  the  range 
3  to  8.  Rased  on  comp orisons  wish  T>ohl  ichting's  [3]  corrected1 
rib-roughness  data,  the  drag  coerfi-ient  function  constant  was 
initially  taken  to  bo  CpQ  «  40.  The  predictions  presented  in  this 
work  are  all  for  fully  developed  flow  between  parallel  plates  or  in 
pipes.  The  motivation  for  using  these  conditions  is  the  great  wealth 
of  experimental  data  that  exists  for  these  configurations.  The 
boundary-layer  equations  (1)  and  (2),  are  readily  reduced  for  these 
conf igurat ions .  All  of  the  comparis^  made  below  are  for  rectangu¬ 
lar  rib  roughness . 

Webb,  Eckert  and  Goldstein  [26]  reported  data  for  fully 
developed  flow  in  square-rib-roughened  pipes  with  a  variety  or  rib 
he ight-to-d iameter  ratios,  k/D,  and  pitch-to-h«ight  ratios,  L/k. 
Figure  6  shows  the  comparison  of  the  predictions  with  their  data  for 
L/k  =  10  and  k/D  =  0.01,  0.02  and  0.04,  jsing  'try-,  =  40.  Inspection 
of  the  figure  reveals  that  the  agreement  is  good  for  all  cases. 

All  of  the  comparisons  above  were  for  rib  ;  itch-to-he  ight.  ratios 
of  L/k  «  10.  Webb  et  al.  also  presented  dot  i  for  L/k  =  20  and  L/k  « 
40.  Figure  7  presents  comparisons  between  predictions  and  the  data. 
At  first  the  same  value  of  Cpn  =  40  was  used  for  the  predictions. 
However,  for  the  mere  widely  spaced  ribs  this  value  of  Cqq  did  not 


^chl icht ing's  [31  rib-roughness  data  was  corrected  in  the  same  manner 
as  that  reported  by  the  authors  [ 4 J  for  his  distributed  roughness. 


provide  satisfactory  results.  As  shown  in  the  figure,  a  value  of  C[>q 
=  20  was  found  to  give  good  results  for  L/k  =  20,  and  a  value  of  Cqq 
=  10  was  found  to  give  good  results  for  L/k  =  40. 

Berger  and  Whitehead  [24]  presented  data  for  fully  developed 
flow  in  square  rib-roughened  pipes  with  a  height-to-diameter  ratio  of 
k/D  =  0.02.  Pitch-to-diameter  ratios  of  10,  7.2,  5  and  3  were 
used.  Figure  8  shows  comparisons  of  predictions  with  these  data  for 
Cdo  =  40.  Inspection  of  the  figure  reveals  that  the  comparison  is 
good  for  the  larger  Reynolds  numbers  where  Cf  is  essentially 
constant.  For  the  more  closely  spaced  roughness,  the  data  show  an 
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Figure  8.  Comparison  with  Data  (denoted  by  symbols)  of  Berger  and 
Whitehead  [24]  for  k/D  =  0.02:  (1)  L/k  =  7.2;  (2)  L/k  = 
10;  (3)  L/k  =  5;  (4)  L/k  =  3;  Cpg  =  40  for  all  Cases 

almost  constant  Cf  for  all  Reynolds  numbers  while  the  predictions 
still  indicate  that  Cf  decreases  with  decreasing  Reynolds  number. 
The  predictions  indicate  the  correct  trends  for  the  effect  of 
pitch-to-he ight  ratio.  The  predicted  value  of  Cf  increases  as  L/k 
increases  from  3  to  7.2  and  then  decreases  slightly  for  L/k  »  10. 


Stukel,  Hopke  and  Nourmohammadi  [25]  presented  data  for  fully 
developed  flow  in  square-r ib-roughened  pipes  for  a  wide  variety  of 
pitch-to-height  ratios  ranging  from  L/k  =  2  to  25.  The  data  shown  in 
Figure  9  are  for  a  height-to~diameter  ratio  of  k/D  =  0.072.  The 
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Figure  9.  Comparison  with  Data  (denoted  by  symbols)  of  Stukel  et 
al.  [25];  Re  =  150,000;  k/D  =  0.072;  CD0  =  10  for  all 
Cases 

value  of  C[)o  used  for  the  predictions  shown  was  Cpg  =  10.  Inspection 
of  the  figure  reveals  that  the  model  predicts  the  correct  trends  of 
the  influence  of  L/k.  However,  the  model  predicts  the  maximum  value 
of  Cf  at  L/k  =  8  while  the  data  indicate  a  maximum  near  L/k  =  12. 

Han,  Glicksman  and  Rohsenow  [12]  presented  data  for  fully 
developed  flow  in  square-rib-roughened  channels  with  a  wide  variety 
of  L/k  and  k/Df,  values.  Comparisons  with  selected  data  are  shown  in 
Figures  10  and  11.  Figure  10  shows  comparisons  for  L/k  =  5  and  k/D^ 

=  0.076  and  O.O1^.  The  value  used  for  C^q  in  Figure  10  was  *10.  From 
the  figure  it  is  seen  that  the  predictions  agree  with  the  data  for 
Reynolds  numbers  above  approximately  10,000  where  both  the  data  and 
predictions  shown  an  essentially  constant  value  of  Cf.  For  lower 
Reynolds  numbers  the  model  predicts  a  more  rapid  decrease  in  Cf  than 
the  data.  Figure  11  shows  comparisons  for  k/D^  =  0.056  with  L/k  =  10 


and  20.  Following  the  experience  with  the  data  of  Webb  et  al.  for 
L/k  =  10,  Cqo  was  taken  to  be  40;  and  for  L/k  =  20,  Cpg  was  taken  to 
be  20.  Inspection  of  the  figure  reveals  that  the  model  overpredicts 
the  skin  friction  coefficient  for  these  data  by  about  20$. 
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Figure  10.  Comparison  with  Data  (denoted  by  symbols)  of  Han  et  al. 

[12]  for  L/k  =  5;  (1)  k/Dh  =  0.076;  (2)  k/Dh  =  0.046;  CD0 
=  40  for  all  Cases 
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Figure  11.  Comparison  with  Data  (denoted  by  symbols)  of  Han  et  al. 

[12]  for  k/D  =  0.056;  (1)  L/k  =  10,  CD0  =  40;  (2)  L/k  = 


SECTION  VI 


DISCUSSION  AND  CONCLUSIONS 

From  the  comparisons  shown  in  the  preceding  section,  it  is  seen 
that  the  model  reproduces  important  salient  features  of  the  turbulent 
flow  over  two-dimensional  rib  roughness.  However,  it  is  found  that 
the  discrete  element  approach  has  a  much  narrower  range  of  applica¬ 
tion  for  different  two-dimensional  surface  geometries  than  the  very 
broad  range  which  has  previously  been  demonstrated  by  the  authors  [1, 
17]  for  three-dimensional,  distributed  roughness.  This  is  shown  by 
the  necessity  of  changing  the  drag  coefficient  parameter,  Cdo»  for 
different  values  of  L/k.  This  should  not  have  been  surprising  since 
the  length  of  attached  stream-wise  flow  influences  the  shape  of  the 
upstream  velocity  profile  (momentum  distribution)  which  determines  in 
part  the  pressure  distributions  on  the  upstream  and  downstream  faces 
of  a  two-dimensional  roughness  element.  The  details  of  the  nature  of 
the  windward-and-leeward  separated  regions  and  the  interaction  with 
the  main  streamwise  flow  are  important  in  establishing  the  form  drag 
of  the  roughness  element. 

Based  on  the  comparisons  with  experimental  data,  the  following 
values  of  Cqq  are  proposed: 

CDo  =  'tO,  L/k  <  10  (12) 

CDo  =  '100/ (L/k),  L/k  >  10 

This  formulation  gives  reasonable  agreement  for  Re  >  10,000  for  three 
of  the  four  data  sets  considered.  The  exception  is  the  data  of 
Stukel  et  al.  [25].  In  this  data  set  the  rib  height  was  1 of  the 
pipe  radius,  and  the  axial  variation  of  the  centerline  velocity  was 
significant  between  the  ribs.  This  effect  was  not  taken  into  account 
in  this  study. 

Although  the  two-dimensional  rib-roughness  model  developed 
herein  was  successful  in  predicting  many  important  salient  features 
of  the  dependence  of  skin  friction  on  rib  spacing,  rib  size,  and 
Reynolds  number,  overall  the  model's  results  were  not  nearly  as 
satisfying  a3  those  of  the  previously  developed  three-dimensional 
roughness  model.  With  the  three-dimensional  distributed  roughness, 


21 


the  averaging  processes  employed  in  the  x-z  and  y-z  planes  are 
physically  more  meaningful  and  appropriate  than  for  the 
two-dimensional  roughness. 

Because  equations  (1)  and  (2)  are  the  correct  equations  for 
discrete  element  considerations,  any  further  advance  in  the 
predictive  ability  of  the  technique  developed  herein  must  come 
through  better  modeling  of  the  blockage  factors  and  form  drag.  The 
two-dimensional  roughness  elements  are  in  discrete  locations, 
perpendicular  to  the  streamwise  flow,  and  cause  significant  local 
perturbations  including  streamline  curvature.  The  key  to  modeling 
two-dimensional  roughness  in  the  discrete  element  averaging  sense  is 
specification  of  the  separation  lengths  and  the  form  drag  coefficient 
variation  over  the  roughness  geometry.  Since  the  discrete  element 
technique  and  two-dimensional  rib-roughness  models  are  attempts  to 
construct  a  completely  parabolic  (marching)  methodology  for 
situations  which  in  reality  possesses  significant  embedded  elliptic 
regions,  one  should  not  expect  a  two-dimensional  boundary  layer 
approach  with  any  roughness  model  to  have  broad  applications  for 
turbulent  flow  over  r ib-roughened  surfaces. 

Another  area  of  importance  for  further  investigation  is 
two-dimensional  r ib-roughness  not  oriented  perpendicular  to  the 
streamwise  flow.  Cross  flow  effects  for  this  condition  are  likely  to 
be  important,  and  the  model  developed  herein  will  need  further 
modification. 
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